Absence of residual quasiparticle conductivity in the underdoped cuprate YBa2Cu408 
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We report here measurements of the in-plane thermal 
conductivity k of the stoichiometric underdoped cuprate 
YBa2Cu408 (Y124) below IK. k{T) is shown to follow a sim- 
ple, phononic T'^ dependence at the lowest T for both current 
directions, with a negligible linear, quasiparticle contribution. 
This observation is in marked contrast with behavior reported 
in optimally doped cuprates, and implies that extended zero- 
energy (or low-energy) quasiparticles are absent in Y124. 



Experimental evidence for rf-wave superconductivity in 
high-Tc cuprates is now well established. The presence 
of nodes in the gap is expected to produce a finite den- 
sity of well-defined quasiparticle (QP) excitations at low 
energies that dominate the low-T physics. For a pure 
d-wave superconductor with line nodes on the Fermi sur- 
face, for example, the density of states (DOS) is linear 
in the excitation energy, giving rise to a dependence 
of the low- T specific heat and thermal conductivity (as- 
suming a constant scattering rate). This excitation spec- 
trum, however, is altered significantly in the presence of 
impurities. For an anisotropic two-dimensional (2D) su- 
perconductor with scattering in the unitary limit, a band 
of impurity states is expected to develop whose width 7 
grows with increasing impurity concentration riimp, lead- 
ing to a finite zero-energy DOS Lee showed that 
the residual conductivity is independent of Uimp, the re- 
sult of compensation between the increased DOS and a 
reduction in the associated transport lifetime This 
residual, or "universal" conductivity develops at low T 
in the so-called "dirty" limit, fcsT < 7. Similarly, the 
low- T k{T) will be dominated by QP states in the vicin- 
ity of the line nodes, and is given by the expression |^ 

KrejT^in/d)ikl/3)ivF/v2) (1) 

where n/d is the stacking density of Cu02 planes and vp 
and V2 are the energy dispersions (QP velocities) perpen- 
dicular and tangential to the Fermi surface respectively. 



The issue of localization of electronic states in rf-wave 
superconductors, however, is a complex problem and 
many competing viewpoints prevail [^-^j. Balatsky and 
Salkola Q, for example, argue that the long-range na- 
ture of hopping between impurity states along the nodal 
directions leads to strong overlap of the impurity wave 
functions along the diagonals of the square lattice, and 
ultimately to "extended" impurity states. Senthil and co- 
workers [^,0, on the other hand, argue that quantum in- 
terference effects destabilize the extended QP states and 
lead to a vanishing DOS at zero energy (See also [H for 
a similar conclusioii). They coined the term "supercon- 
ducting insulator" |7| to describe such a superconductor 
with localized states. 

Experimentally, the observation of a finite linear term 
in k{T) in pure and Zn-doped YBa2Cu307_5 (Y123) by 
Taillefer et al. appeared to confirm the existence of 
zero-energy quasiparticles. Moreover, the size of this 
term was indeed found to be "universal", i.e. inde- 
pendent of Zn concentration, in agreement with Lee's 
prediction. Similar behavior was later reported for 
Bi2Sr2CaCu208 (Bi2212) by Behnia et al., using irra- 
diated crystals More recently, the magnitude of 
Kres/T in Bi2212 (« 0.15 mW/cm.K^) was shown to be 
consistent with absolute values of vf/v2 estimated from 
angle-resolved photoemission (ARPES) pl] , p^ . 

Despite this apparent consistency between theory and 
experiment, Kres/T has only been reported for two com- 
pounds, both at their optimum doping level, and it is 
not immediately obvious how Kres/T will vary across the 
phase diagram. ARPES and penetration depth mea- 
surements [ p^ support claims that vf/v2, and thereby 
Kres/T, increase in the underdoped regime. On the other 
hand, in certain underdoped cuprates, where Te has been 
suppressed in high magnetic fields 0,|5|, there is a 
marked tendency towards localization below Te, suggest- 
ing a vanishing QP contribution at low T. Clearly, low- T 
k(T) measurements on underdoped cuprates are impor- 
tant to help clarify this seemingly contradictory behavior. 



1 



With this in mind, we have carried out the first low- 
T k{T) measurements on the underdoped cuprate Y124 
(Tc = 80K), which is a self-doped, stoichiometric cuprate 
and therefore relatively free of disorder. Below 0.25K, 
k(T) « for both a- and &-axis currents, consistent 
with a phonon heat conduction in the ballistic regime. 
The residual linear QP term, however, is either absent 
or is negligibly small, with an upper bound of 0.02 
mW/cm.K^. This result reveals that the universal con- 
ductivity scenario breaks down dramatically in under- 
doped Y124. One compelling possibility is that the QP 
states in Y124 are localized at low T, due to the proxim- 
ity to the superconductor/insulator (S/I) boundary, and 
therefore do not contribute to the low-T heat transport. 

The Y124 crystals were grown by a flux method de- 
scribed elsewhere For this particular study, three 
plate-like crystals were selected, two with their longest 
dimension along the &-axis, the other along the a-axis. 
Approximate dimensions were 0.25 x 0.16 x 0.015 mm^ 
for the Ka crystal (labelled hereafter as ajjl) and 1 x 0.09 
X 0.05 mm^ and 0.8 x 0.25 x 0.06 mm^ for the two Kb 
crystals, b(jl and bp. Tc = 80K for all crystals, with a 
transition width, measured resistively, of less than IK. 

k{T) for each crystal was measured between 0.14K and 
IK using a conventional steady-state four-probe tech- 
nique that allowed the electrical resistivity Pa.b{T) of each 
sample to be measured in situ without changing the con- 
tact configuration. Gold wires were attached as electrical 
contacts using Dupont 6838 silver paint. The Pa,b{T) be- 
havior was found to be in excellent agreement with pre- 
vious measurements [ p^ , with room temperature values. 
Pa = 350 /ificm and pt, = 90 pQ.cva (for both 6-axis crys- 
tals). This large in-plane anisotropy arises from the high 
conductivity of carriers on the quasi-lD CuO chains that 
run parallel to the &-axis (see schematic inset to Figure 
1) and confirms not only the high quality of the crystals 
used in this study, but also that current flow in each case 
is uniaxial. The temperature gradient was measured by 
two Ru02 thermometers connected to the "voltage" con- 
tacts through the gold wires, and the thermometers were 
supported by long, thin superconducting Nb-Ti wires to 
minimize heat losses. Uncertainties in the absolute mag- 
nitudes of k(T) (p(T)), mainly due to the finite contact 
dimensions on these small crystals, are estimated to be 
around 15% for Kf, (pb) and around 25% for Kq {pa)- 

The n{T) data for all crystals are shown on double- 
logarithmic axes in Figure 1. The variation of K,h{T) 
for the two fe-axis crystals is almost identical over the 
whole temperature range studied, giving us confidence in 
the reproducibility of our data. Below 0.25K, Ha and Kb 
both vary approximately as T^, consistent with phonon 
heat transport in the boundary-scattering limit. Above 
0.25K, Ka{T) deviates more strongly from a depen- 
dence. The origin of the enhancement of kjj over Ka above 
T = 0.25K is not understood at present, though we as- 
sume it reflects an additional contribution to k from the 



CuO chains; either QP conductivity on the chains devel- 
ops swiftly above 0.25K (note that the CuO chains, being 
quasi-lD, may be susceptible to charge ordering at very 
low T"), or there exists an additional channel for phonon 
heat propagation jl^ along the chains that reduces the 
effects of phonon scattering beyond the ballistic regime. 
Further measurements in a magnetic field are envisaged 
to clarify the origin of this anisotropy. 

In order to look for evidence of a linear K^es, we have 
re-plotted the low T data in Figure 2 as Ka,b/T versus 
and fitted each data set below 0.25K to the expression 
Ka,b ^AT + BT^ |l|]. The coefficients for each fit are A 
= 0.011, - 0.007 and 0.006 (± 0.02) mW/cm.K^ and B 
= 6.50, 7.00 and 10.67 (± 0.50) mW/cm.K^ for ajjl, bjjl 
and bp respectively. 

In the boundary-scattering limit, Kph is given by 

Kph = 1/3/3 < vph > loT"" (2) 

where /3 is the phonon specific heat coefficient, < Vph > 
the average acoustic sound velocity and Iq = 2w/y/TT is 
the maximum phonon mean free path. Here, w represents 
a mean width of the rectangular-shaped crystal. Taking 
the dimensions of our crystals and suitable values for /3 
(= 0.5 ± 0.1 mJ/mol.K^) [|o) and < Vph > 5 ± 1 
X 10^ cm/s) ||2|], we obtain estimates for Kph/T^ — 4.1 
± 1.2, 5.25 ± 1.5 and 9.55 ± 2.0 mW/cm.K^ for ajjl, 
bjjl and bp respectively Given the uncertainties in 
measuring dimensions and contact distances, we believe 
these values compare favourably with the experimental 
values. More importantly, the size of the term for the 
two Z)-axis crystals scales well with w and we conclude 
that the contribution is indeed simply the phonon 
contribution in the ballistic regime. 

The most striking result here is the complete absence 
(to within our experimental accuracy) of the residual lin- 
ear term in the low-T i^{T) for both chain and plane cur- 
rent directions. It should be emphasized, of course, that 
a zero linear term in K}, also implies a negligible Kres 
within the planes, meaning we have effectively confirmed 
the absence of the universal QP term in Y124 in all three 
samples. Moreover, for there to be any finite zero- T in- 
tercept in k/T, it would require Kph{T) below 0.15K to 
vary as T'''+" with n > 0, which is simply not physical, 
given that the lattice heat capacity is strictly cubic below 
IK. Thus, we are confident that the main result of this 
Letter, namely the absence of Kres in Y124, is robust. 
For comparison, we also show in the inset to Figure 2, 
Kab/T for optimally doped Bi2212 [lO| ] measured with the 
same experimental set-up p^ . In Bi2212, we can clearly 
distinguish a finite Kres/T ~ 0.15 mW/cm.K^ (shown by 
a dotted line) , that is an order of magnitude larger than 
the upper limits for Kres/T in Y124. 

Despite the overwhelming case for d-wave pairing in 
high- Tc cuprates, there is still limited, direct evidence 
for a d,j.2_y2 order parameter in Y124. Thus, before dis- 
cussing our result in terms of nodal QP states, we should 
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first examine the possibility that there is a finite gap ev- 
erywhere on the Fermi surface in Y124. First of all, the 
orthorhombic distortion in Y124, induced by the chains, 
introduces some d + s admixture in the gap function. 
As the s-component is increased from zero, the posi- 
tion of the nodal lines are first shifted away from (7r,7r), 
but as the s-component becomes comparable with the 
d-component, a nodeless gap may form. Secondly, mag- 
netic impurities are thought to induce a local imaginary 
component that could also give rise to a fully gapped 
state and a suppression of low T thermal transport . 
This latter possibility, however, is not supported by spe- 
cific heat data taken on crystals from similar batches to 
those studied here, which show no sign of a low T Schot- 
tky anomaly arising from such magnetic impurities [^ . 
In addition, power law penetration depths have now been 
observed down to 2K for both the a- and 6- axes in Y124 
|^,|6|, suggesting a simple nodal gap picture is equally 
applicable to Y124. Of course, we cannot rule out com- 
pletely the possibility of a finite gap in Y124, but if it does 
exist, it would have to be vanishingiy small. In what fol- 
lows, therefore, we assume that nodal lines are present in 
Y124 and turn to consider what might be happening to 
the low-energy QP states in their vicinity. 

From (1), we see that Kres/T is directly proportional 
to the ratio vf/v2 at the nodal positions, so a negligible 
Kres/T may indicate a sharp gap feature at the nodes at 
very low energies, induced either by doping, impurities 
or structural modifications. However, as noted above, 
vfIv2 is expected to increase as we move into the un- 
derdoped regime, and indeed, independent estimates of 
vp/v2 from the slope of the low-T penetration depth p^ ] 
in Y124 yield an estimate for Kres/T that is larger than 
those for both optimally doped Bi2212 and Y123. More- 
over, given that band structure estimates for vp are sim- 
ilar for Y123 and Y124 even with our upper bound 
estimate for Kres/T (= 0.02 mW/cm.K^), we obtain a 
physically unrealistic value (vp /v^ ~ 2) for the gap slope 
within the nodes in Y124. It appears unlikely, therefore, 
that the gap structure can itself explain a value of Kres/T 
one order of magnitude lower than in Y123. 

Another important consideration is the size of the im- 
purity band, 7. We recall that in the unitary limit, the 
universal conductivity regime develops below /cgT < 7. 
Thus, the observation of Kres depends not only on the 
temperature range of the experiment, but also on the 
energy scale of the impurity band, imposed by the scat- 
tering phase shift. However, taking values for the 6-axis 
residual resistivity po = 0-5 [iSXcm jl^ and plasma fre- 
quency Wp = 2.5eV p^ , we obtain a lower bound es- 
timate of 7 in the unitary limit of 14K, some two or- 
ders higher than the base temperature of our measure- 
ments. (Unfortunately, similar analysis for J/ /a cannot 
be performed due to the difhculties in estimating po from 
Pa{T)). Even in the opposite (Born) limit, where 7 be- 
comes exponentially small and therefore, may fall below 



our measurement range, the product of the DOS and 
the lifetime is also energy independent above ksT > 7. 
Hence, Kres/T should still be constant and finite and the 
same arguments still apply. 

Given these simple yet rather compelling arguments, 
we are left to consider how QP localization might ac- 
count for the absence of Kres/T in underdoped Y124. 
As mentioned above, high magnetic field measurements 
on the underdoped cuprates La2-a;Sra;Cu04 and 
Pr2_a;Cej;Cu04 revealed that pabiT), though " metal- 
lic" at high T, tends towards localization as T OK. 
The origin of this localization phenomenon is unknown 
at present. However, if the field- induced destruction of 
superconductivity leads directly to an insulating phase, 
then it is not unreasonable to assume that this transi- 
tion is from a superconductor with already localized QP 
excitations [Q. The crossover from metallic to insulat- 
ing behavior in the normal state (i.e. above He2) as we 
approach the Mott insulator, suggests an increasing role 
of long-range interactions on the mobility of low-energy 
quasiparticles. The superfluid condensate, suppressed in 
the vicinity of an impurity, becomes ineffective in screen- 
ing completely the Coulomb repulsion between quasipar- 
ticles in the bound state |^ . As we approach the parent 
insulator, we expect such interactions to grow, leading to 
insulating behavior of the quasiparticles above the super- 
conducting condensate and a negligible QP contribution 
to the low-T heat transport |Q. Such localization in a 
nominally clean superconductor is an exciting prospect, 
and measurements on Zn-dopcd or irradiated Y124 are 
envisaged to investigate this possibility further. We note 
here that most impurity models for d-wave (cuprate) su- 
perconductivity fail to take into account the developing 
role of long-range interactions between quasiparticles as 
the S/I boundary is approached. We hope, therefore, 
that this work stimulates renewed theoretical efforts to 
understand the nature of QP excitations in the Cu02 
planes, deep inside the superconducting state, and in par- 
ticular in stoichiometric crystals on the underdoped side 
of the phase diagram. 

In conclusion, we have measured the low- T k{T) of sto- 
ichiometric, underdoped Y124 and have found that, in 
marked contrast to optimally doped Y123 and Bi2212, 
the "universal" QP conductivity term is absent. We 
have considered several interpretations of this intriguing 
result, including localization of the quasiparticles them- 
selves, due to enhanced long-range interactions as we ap- 
proach the S/I boundary. Prior to this work, the obser- 
vation of the universal conductivity in Y123 and Bi2212 
had been widely regarded as solid support for the pic- 
ture of long-lived quasiparticles above the superconduct- 
ing ground state of high-Tc cuprates. Our surprising 
result offers important and timely counter evidence that 
the residual conductivity term is non-universal, and we 
hope it encourages further debate and investigation into 
this critical, and still controversial, issue. 
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Figure Captures 

Fig.l. Ka{T) and ^^(r) of Y124, plotted on double- 
logarithmic axes. The dashed line represents the de- 
pendence expected for phonon heat transport in the bal- 
listic regime. The inset shows a simple schematic of the 
crystal structure of Y124. 

Fig. 2. k/T versus below 0.4K for ajjl (closed cir- 
cles), bjjl (open circles) and bp (open squares). Fits 
to the expression k = AT + BT^ below 0.25K are indi- 
cated by dashed lines for bjjl and bp and by a solid line 
for ajjl. The dotted line represents the universal con- 
ductivity limit for Bi2212. Inset: Comparison of K,a/T 
for Y124 (closed circles) and Kab/T of Bi2212 Q (open 
diamonds), measured in the same apparatus. 
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